U.S.  Army  Workshop  at  Imperial  College: 


IMPERIAL  COLLEGE 

Mechanical  Engineering  -  Room  703 

Introductory  Remarks 


Thursday.  26  June.  1990 

9.00  Introductory  Remarks 
Unsteady  and  two-phase  flows 


R.  Reichenbach 

Chairman 

9.15 

9.45 

A  Brown  K 

M.  Sommerfeld 

Unsteady  fluid  dynamics 

Shock  wave  propagation  in  a  dusty  gas. 

10.15 

Coffee  break 

10.45 

11.15 

11.45 

12.15 

A.D.  Gosman 

A.M.P.K.  Taylor 

H.  McDonald 

M.  Oldfield 

A  calculation  method  for  unsteady,  two-phase  flows. 
Methods  for  measuring  two-phase  flow  properties 
Calculation  of  two-phase  flows. 

Unsteady  pressure  and  heat-transfer  measurements. 

13.00 

Lunch 

Gun  flows 

A.  Brown 

Chairman 

14.30 

15.00 

15.30 

G.  Keller 

R.  Heiser 

Coflee  break 

Interior  ballistic  effects  of  kinetics,  grain  fracture  and 
viscosity. 

Interior  ballistic  flow  modelling  at  the  Ernst  Mach 
Institute. 

16.00 

16.30 

17.00 

G.  Klingenberg 

R.  Perez-Ortiz 

B.  Lawton 

Unsteady  flow  and  combustion  phenomena. 

Velocity,  pressure  and  heat  transfer  in  a  gas  gun 
Instantaneous  heat  transfer  in  gun  barrels. 

19.00  Dinner 


Friday.  29  June  1990 


Unsteady  and  two-phase  flows 


G.  Bergc'.as 

Chairman 

9.00 

9.30 

10.00 

K. C.  Schadow  and  E.  Gutmark 
M.N.R.  Nina 

L.  De  Luca 

Passive  and  active  combustion  control. 

Passive  control  of  combustion  oscillations. 

Burning  stability  of  solid  propellants. 

10.30 

Coffee  break 

11.00 

11.30 

12.00 

12.30 

F.E.C.  Culick 

E.J.  Mularz 

J.A.G.  Aston  and  A.D.  Crowley 
T.  Cebeci 

Pressure  oscillations  in  liquid-fuelled  propulsion 
systems. 

Report  oOn  status  of  hydrogen  fuel-shear-layer 
experiment. 

On  the  flow  structure  in  two-phase  shock  lubes. 
Calculation  of  unsteady  flows  with  separation. 

13.00 

Lunch 

Unsteady  and  two-Dhase  flows 

14.30 

15.00 

C.  Zoltani 

L.  Carr 

Investigation  of  confined,  two-phase,  subsonic  jet  flow. 
Experiments  with  oscillating  airfoils. 

Gun  Hows 

R.  Ogorkiewicz 

Chairman 

15.30 

16.00 

F.  Seiler  and  K.  Zimmerman 

M.  Mach 

A  study  on  gun  ercclo.i  in  a  2Cmm-caiibe>  g^r.  ;„be. 
Studies  on  secondary  combustion  of  gun  exhaust  plumes. 

16.30 

Coffee  Break  and  Discussion. 

17.30 

Closure 

1 


U.S.  Army  Workshop  at  Imperial  College,  28-29  June  1990 


Unsteady  Fluid  Dynamics  Research  at  RMCS  Shrivenham 

A  Brown*,  M  V  Finnisb,  D  Bray°,  K  Knowlesc  ,  G  M  Moss1*  and  S  B  Ellis* 

School  of  Mechanical,  Materials  and  Civil  Engineering 
Royal  Military  College  of  Science,  Shrivenham,  Swindon,  SN6  8LA 


Abstract 

Unsteady  fluid  dynamics  and  aerodynamics  is  an  area  of  major  interest  at  RMCS. 
Within  this  general  area  research  covers  boundary  layer  transition  on  flat  and 
curved  surfaces  in  the  presence  of  pressure  gradients,  jet  aerodynamics,  vortex 
aerodynamics  and  ballistics.  All  of  this  work  is  relevant  to  the  performance  of 
military  equipment,  particularly  gas  turbines,  helicopters  and  jet-lift  aircraft, 
guided  weapons  and  tube-launched  projectiles. 

In  laminar  boundary  layers  it  is  generally  accepted  that  under  zero-pressure 
gradient  conditions  instabilities  can  occur  at  a  critical  point  where  the 
Reynolds  number  based  on  momentum  thickness  has  a  value  of  about  160.  As  this 
Reynolds  number  reaches  200,  turoulent  spots  may  occur  which  grow  and  lead  to  a 
fully  turbulent  boundary  layer.  Freestream  turbulence  plays  a  part  in  the  onset 
and  extent  of  transition,  as  does  surface  vibration.  Recent  work  has  investigated 
the  relationship  between  the  frequency  and  amplitude  of  surface  vibrations  and 
the  onset  of  transition1.  Of  particular  interest  in  turbomachinery  applications, 
however,  is  the  influence  of  pressure  gradients  on  boundary  layer  transition. 
Extensive  experimental  and  theoretical  work  has  been  carried  out  for  both  adverse 
and  favourable  pressure  gradients.  The  work  was  aimed  at  a  better  understanding 
of  heat  transfer  and  skin  friction  on  turbine  blade  surfaces  but  is  equally 
applicable  to  flow  over  any  curved  surface. 

In  the  case  of  convex  curved  surfaces,  the  most  successful  theoretical  approach 
has  been  found  to  be  a  modification  of  the  STANS  program  from  Crawford  and  Kays'2. 
An  important  characteristic  of  any  boundary  layer  prediction  program  is  the 
transition  function:  the  original  version  of  STAN5  contained  a  transition 
function  of  a  sinusoidal  nature.  This  could  be  called  into  action  at  a  pre¬ 
specified  critical  momentum- thickness  Reynolds  number,  usually  at  a  value  of  200. 
At  this  point  the  turbulent  stresses  were  gradually  increased  from  zero  to  their 
fully-turbulent  value  at  the  point  where  the  momentum-thickness  Reynolds  number 
was  equal  to  twice  the  critical  value .  Roberts  and  Brown3  developed  an 
alternative  transition  function  based  on  a  velocity  gradient  factor4-5  and  the 
total  strain  rate  up  to  the  position  of  start  of  transition.  This  fitted  well 
with  their  measurements  on  the  suction  surface  of  a  gas  turbine  blade  in  a 
cascade.  The  foregoing,  when  coupled  with  a  function  of  freestream  turbulence 
intensity,  helped  in  predicting  heat  transfer  rates  between  the  hot  gases  and  the 
fade  surface. 


a  Professor,  Chairman  of  School  and  Head  of  Thermal  Power  Group, 
b  Research  Officer,  Thermal  Power  Group, 
c  Lecturer,  Aeromechanical  Systems  Group, 
d  Senior  t»cfv-er.  Aercccoh*r.i''a!1  S’/rte’ —  Group, 
e  Senior  Research  Officer,  Aeromechanical  Systems  Group. 


< 

J 


2 


On  concave  surfaces  the  boundary  layer  is  subjected  to  an  instability  due  to 
centrifugal  forces.  This  results  in  a  secondary  flow  in  the  form  of  streamwise 
counter-rotating  vortices.  These  Goertler  vortices  are  known  to  be  a  modulator  of 
transition  to  turbulence  and  the  boundary  layer  prediction  programs  currently 
available  have  a  poor  performance  in  such  flows.  Experimental  and  theoretical 
work  is  currently  in  hand  to  study  this  flow  pattern0-7.  Predictions  of 
streamwise  velocity  profiles  and  ranges  of  vortex  amplitude  agree  well  with 
experimental  measurements  for  zero-pressure  gradient.  Predictions  for  Falkner- 
Skan  flows  indicate  an  increased  stability  to  Goertler  vortices  for  favourable 
pressure  gradients  and  decreased  stability  for  adverse  gradients. 

Jet  aerodynamics  research  at  RMCS  is  currently  concerned  with  two  topics: 
propulsion  nozzle  mass  flow  control  using  swirl  and  impinging  jets  in  cross-flow. 
The  former  has  involved  an  analytical  study®  which  has  demonstrated  the 
feasibility  of  the  scheme0;  this  is  not,  however,  an  unsteady  flow  phenomenon. 
Impinging  jets,  on  the  other  hand,  are  a  well-known  source  of  flow 
unsteadiness10. When  a  jet  impinges  on  a  surface  normal  to  its  axis  a  wall  jet  is 
formed  which  spreads  out  radially  from  the  impingement  point.  If  a  cross-flow 
(parallel  to  the  wall)  exists  then  the  wall  jet  will  eventually  stagnate  and  roll 
up  to  form  a  vortex,  arranged  as  a  horse-shoe  about  the  jet  (see  Fig  1).  This 
flow-field  is  of  particular  interest  in  the  development  of  V/STOL  aircraft,  where 
the  cross-flow  may  be  due  to  ambient  wind  or  aircraft  motion.  With  multiple  jets 
there  is  the  added  complexity  of  the  fountain  formed  by  the  two  opposing  wall 
jets. 

The  mean  position  and  strength  of  the  ground  vortex  are  important  due  to  their 
effects  on  aircraft  pressure  loads  and  air  intake  flows  (where  there  is  the 
specific  problem  of  hot  gas  ingestion).  Some  of  the  factors  which  affect  the 
ground  vortex  position  have  been  investigated  previously  but  there  is 
considerable  scatter  in  the  experimental  data  and  very  little  high  jet  Mach 
number  data.  The  current  project 11  has  been  looking  at  the  influence  of  nozzle 
height  (2  to  S  diameters),  nozzle  pressure  ratio  (1.05  to  4),  cross-flow- to-jet 
velocity  ratio  and  vector  angle  (+10  to  -20  degs)  for  single-  and  twin-jet 
configurations  with  both  fixed  and  moving  ground  planes  (simulating  ambient  wind 
or  aircraft  motion).  For  the  twin-jet  configurations  the  effect  of  nozzle  splay 
is  also  being  investigated.  The  programme  of  work  is  both  experimental  and 
numerical,  a  particular  aim  being  to  determine  the  difference  in  ground  vortex 
position  with  a  moving,  rather  than  fixed  ground  plane.  This  information  is 
needed  before  wind  tunnel  tests  are  conducted  on  specific  V/STOL  aircraft 
configurations  in  ground  effect. 

An  extensive  series  of  experiments  is  being  conducted  in  the  RMCS  low-speed, 
open-jet  wind  tunnel,  which  is  equipped  with  a  rolling-road  (1.6m  x  lm)  and  a 
compressed  air  supply  (to  single  or  twin  25mm  diameter  nozzles).  Tests  have 
consisted  of  pitot-static  probe  traverses  a  few  millimetres  above  the  ground 
plane,  upstream  of  the  ground  vortex  in  the  plane  of  symmetry.  This  has  allowed 
the  vortex  penetration  to  be  quantified  (see  Fig. 2)  and  plotted  against 
parameters  such  as  cross-flow- to-jet  velocity  ratio  (Fig. 3).  Comparisons  with 
fixed  ground  board  statics  showed  this  approach  to  be  perfectly  satisfactory 
(except  directly  under  the  vortex  core  where  the  probe,  not  surprisingly,  gave  a 
lower  static  pressure).  A  correlation  has  been  derived  for  the  positions  of 
ground  vortex  core,  separation  point  and  maximum  penetration  point12  (see  Fig. 2). 

Flow  visualisation  has  been  conducted  using  smoke  and  laser  light-sheet.  Results 
for  a  twin-nozzle  arrangement  have  shown  the  existence  of  a  "spike"  in  the  plane 
of  symmetry.  The  ground  vortex  formed  with  twin  nozzles  is  seen  to  be  even  more 
unsteady  than  that  due  to  a  single  nozzle  and  the  flow-field  seems  to  be  very 
sensitive  to  cross-flow  yaw.  Computational  work  to  date13  has  used  the  PHOENICS 
finite-volume  code  to  perform  2-D  calculations  with  the  standard  k-e  turbulence 
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model.  Initial  results  are  encouraging  (see  Fig.  4)  and  work  is  proceeding  to 
extend  the  computations  to  the  full  3-D  case. 

Vortex  aerodynamics  has  been  a  subject  of  research  at  RMCS  for  many  years  due 
particularly  to  its  importance  in  guided  weapon  flight  and  in  helicopter  rotor 
performance.  Current  work  is  investigating  rotor/fuseiage ,  blade/vortex  (BVI)  and 
rotor/tail-rotor  interactions.  A  transverse  vortex  generator  has  been  built.  This 
consists  of  a  symmetrical  aerofoil  section  (NACA  0013)  spanning  the  8ft,  low- 
speed  test  section  of  the  RhCS  closed  wind-tunnel  and  capable  of  being  pitched 
rapidly  to  a  modest  angle  of  incidence  (approx.  10°).  This  sheds  a  stc_ting 
vortex  which  then  moves  downstream  over  a  test  body.  The  test  body,  which  could 
be  an  aerofoil  for  BVI  experiments^  or  a  body  of  revolution  for  rotor/fuselage 
interaction  experiments,  is  instrumented  with  high-frequency  pressure  switches 
allowing  pressure  vs.  time  maps  to  be  drawn  for  the  pressure-tapped  area.  The 
interactions  being  modelled  in  these  experiments  are  idealised  events,  not  of 
themselves  representative  of  real  interaction  geometries,  but  capable  of 
providing  data  for  computer  code  validation  or  empirical  input. 

Within  the  area  of  aero-ballistics  RMCS  is  investigating  the  flight  of  tumbling 
projectiles  as  part  of  a  major  study  which  is  intended  to  put  onto  a  proper 
scientific  basis  the  estimation  of  safety  traces  for  firing  ranges.  The  areas 
which  are  currently  being  examined  include  the  following. 

i)  The  ricochet  of  typical  projectile  shapes  from  different  media,  including 
wet  and  dry  sand,  concrete,  pebbles,  turf  and  water.  There  are  several 
aspects  to  these  investigations  which  are  designed  to  highlight  the  physical 
mechanisms  involved  in  the  ricochet  process  and  to  provide  post-ricochet 
data  for  subsequent  trajectory  prediction. 

ii)  The  generation  of  a  full  six-degree-of-freedom  trajectory  simulation  for  a 
complete  range  of  attitudes  and  with  no  small-angle  or  linearising 
assumptions. 

iii)  The  production  of  aerodynamic  data  for  the  above  including  the  effects  of 
round  asymmetry  caused  by  impact  and  aerodynamic  parameters  for  rounds  with 
large  angular  velocities. 

iv)  The  trim  states  for  asymmetric  rounds  having  varying  degrees  of  axial  spin. 

The  attached  figures  give  some  specimen  results.  Figure  5  is  the  predicted  drag 
coefficient  for  an  M80  round  at  a  flight  Mach  number  of  2.5,  whilst  Figure  6 
gives  some  measured  speed  vs  time  data  for  7.62mm  tracer  rounds  after  ricochet. 
The  rate  of  decay  of  these  is  proportional  to  the  effective  drag  coefficient  of 
the  round.  The  inference  is  that,  in  some  cases,  a  tumbling  round  re-stabilises 
into  low-yaw  angle  flight. 

Figures  7  and  8  show  the  predicted  flight  paths  for  a  bent  round  flying  in  a 
dense  medium.  The  changes  in  the  sign  of  the  heading  angle  are  particularly 
noteworthy  -  if  this  were  considered  to  be  the  flight  path  in  the  sand  of  a  stop 
butt,  the  round  could  subsequently  emerge  flying  either  to  the  right  or  to  the 
left  of  the  original  aaimuth  heading. 
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The  shock  propagation  in  a  gas-particle  mixture  exhibit  certain 
characteristics  due  to  the  relaxation  phenomena  associated  wtth 
the  dispersed  particles.  This  results  ir.  a  relaxation  region  be¬ 
hind  the  shock  wave  and  a  decay  of  the  shock  wave  propagation 
into  a  dusty  gas.  The  decay  rate  is  dependent  on  the  shock 
strength  and  the  properties  of  the  particles  e.  g.  i,,2e,  material 
density  and  particle  loading.  Experimental  resuits  showed  that 
the  shock  strength  initially  is  rapidly  decaying  and  then  approa¬ 
ching  an  equilibrium  value.  This  equilibrium  value  iS  only  depen¬ 
dent  on  the  shock  strength  and  the  particle  loading. 

Experiments  and  one-d imer.SiOna  1  numerical  simulations  wiii  be 
presented  which  show  the  above  described  phenomena.  Furthermore, 
results  will  be  discussed  which  consider  two-dimensional  effects 
of  shock  wave  propagation  in  dusty  gases,  e.  g.  boundary  layers 
shock  wave  reflection  at  wedges  and  area  constrictions. 
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ABSTRACT 

An  outline  is  provided  of  CFD  methodology  for  predicting  unsteady,  disperse,  two- 
phase  flows  with  accompanying  heat  and  mass  transfer;  and  sample  applications 
are  presented.  The  method  is  based  on  a  continuum  description  of  the  continuous 
phase  and  a  Lagrangian  model  of  the  discrete  phase,  which  may  consist  of 
particles,  droplets  or  bubbles.  Models  are  included  to  represent  turbulence  effects 
on  both  phases  and  the  processes  of  breakup,  collision  and  coalescence  of  the 
discrete  phase. 

The  governing  equations  are  solved  by  a  fully-i mplicit  algorithm  based  on  the 
existing  PISO  method  for  single-phase  flows.  The  two-phase  version  contains 
provision  for  the  strong  mass,  momentum  and  thermal  energy  intercharges  which 
may  take  place  between  the  phases,  and  also  incorporates  a  temporal  subcycling 
procedure  to  improve  accuracy  and  efficiency  in  circumstances  where  the  discrete 
phase  velocities  are  appreciably  higher  than  those  of  the  continuous  phase. 

The  performance  of  the  methodology  will  be  demonstrated  through  extracts  from 
validation  studies  for  pulsed  liquid  sprays  produced  by  a  pressure  atomiser,  as 
used  m  reciprocating  engines. 
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Abstract 


The  presentation  will  provide  an  overview  of  the  optical  methods 
which  exist  for  the  measurement  of  particle  size,  velocity,  size-velocity 
crcss-correlation  and  particle  flux  and  concentration  in  dispersed 
two-phase  flows.  The  available  instruments  can  be  conveniently 
categorized  into  two  types:  'integral  techniques'  and  'single  particle 
counters'.  A  well-known  example  of  the  integral  technique  is  the 
Fraunhofer-  Effraction  instrument  for  the  measurement  of  particle  size 
distribution  integrated  along  the  direction  of  propagation  of  a  laser  beam. 
The  technique  can  size  in  the  range  between  about  1  and  several  hundred 
pm  at  maximum  concentrations  of  the  dispersed  phase  limited  to 
obscuration  of  about  50%  without  corrections  being  required.  The  single 
panicle  counters  (SPCs)  make  measurements  of  the  size  of  individual 
particles  in  a  test  volume  wnich  is  generally  smaller  than  1  mm^.  SPCs 
are  frequently  combined  with  laser-Doppler  velocimeters  to  provide 
simultaneous  measurements  of  size  and  velocity  and  hence  cross 
correlation  also.  SPCs  can  be  further  sub-divided  into  methods  for 
non-spherical  and  spherical  particles.  Ncn-spherical  particles  can  be 
sized  using  either  the  amplitude  or  the  visibility  of  the  Doopper  signal  in 
the  range  between  about  1  -  several  hundered  pm  and  about  5-100  pm 
respectively.  The  accuracy  of  measurement  depends  on  both  the 
concentration  and  shape  of  the  particles.  Spherical  particles  can  be 
conveniently  sized  by  a  third  variation  of  the  SPC,  namely  a  phase-Doppler 
anemometer.  The  basis  of  this  technique  is  that  the  droplet  acts  as  a 
'magnifying  glass'  for  the  fringes  in  the  test  volume  of  the  anemometer 
and  the  instrument  measures  the  magnification  and  hence  infers  the 
droplet  diametfc  r.  Sizes  in  the  range  from  a  few  to  several  hundered  pm 
can  be  measured  and  at  concentrations  which  have  resulted  in  obscuration 
of  the  beams  in  excess  of  30%.  Some  examples  of  measurements  wi'i  be 
provided. 
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Computational  techniques  used  in  simulation  of  two-phase  flows  can  be  broadly 
categorized  into  two  approaches,  viz.  the  Eulerian-Eulerian  analysis  and  the 
Eulerian-Lagrangian  analysis.  Both  techniques  involve  computing  the  continuous  phase 
using  an  Eulerian  analysis.  The  influence  of  the  particulate  phase  (either  particles  or 
droplets)  on  the  continuous  phase  is  accounted  for  by  inclusion  of  inter-phase  coupling  terms 
in  the  Eulerian  equations,  which  in  the  absence  of  these  terms  would  be  the  usual 
Navier-Stokes  equations.  The  particulate  phase,  on  the  other  hand,  may  be  treated  with 
either  a  continuum  or  a  discrete  model.  The  Eulerian-Eulerian  technique  uses  a  continuum 
model  for  the  particulate  phase  and  is  commonly  termed  the  two-fluid  model.  This  approach 
models  a  dense  granular  bed  very  conveniently  and  this  undoubtedly  accounts  for  its 
popularity  in  modeling  gun  interior  ballistics  where  large  particle  loading  ratios  occur  over 
most  of  the  cycle  (see  for  example  Gough[l]  and  Gibeling  et  al.[2]).  The  second  approach 
employs  Lagrangian  description  to  analyze  the  motion  of  the  particulate  phase  using 
computational  "particles"  to  represent  a  collection  of  physical  particles  (or  droplets). 
Newton’s  law  of  motion  is  employed  to  simulate  the  particle  motion  under  the  influence  of 
the  local  environment  produced  by  the  continuous  phase. 

In  simulation  of  evaporating  sprays,  it  becomes  necessary  to  account  for  the  fact  that 
the  discrete  is  phase  not  mono-dispersed.  To  accomplish  this  in  the  Eulerian-Eulerian 
methodology,  the  two-fluid  model  can  be  generalized  into  a  multi-fluid  model.  However,  the 
CPU  time  requirements  increase  rapidly  with  the  increasing  number  of  particle  size  classes 
since  an  extra  "fluid"  has  to  be  added  for  every  particle  size  class.  Also,  numerical  dissipation 
and  dispersion  could  be  a  major  problem  where  it  might  be  necessary  to  maintain  sharp 
inter-phase  boundaries.  The  Eulerian-Lagrangian  analysis  treats  the  particl  size  as  one  of 
the  attributes  assigned  to  computational  particles  and  hence  has  no  trouble  simulating  flows 
which  involve  changing  particle  size.  Further,  this  approach  involves  integration  of  ODE’s 


for  the  particulate  phase  and  hence  is  numerically  efficieni,  in  addition  to  eliminating  the 
numerical  dispersion.  The  deterministic  nature  of  the  particle  dynamics  allows  incorporation 
of  models  for  turbulent  dispersion,  agglomeration,  collision,  etc. 

In  Eulerian-Lagrangian  algorithms,  the  coupling  terms  for  the  continuous  phase 
equations  can  be  computed  from  the  knowledge  of  the  trajectories  for  the  representative 
particles  and  their  attributes  at  the  intersection  of  the  trajectories  with  Eulerian  cell 
boundaries.  The  coupling  terms  can  also  be  computed  from  the  instantaneous  distribution  of 
the  computational  particles  in  the  domain.  The  former  approach  has  been  employed,  for 
example,  by  Crowe,  et  al.  [3]  and  Gosman  and  Ioannides  [4],  while  the  later  approach  has 
been  used  by  Duckowicz  [5]  and  Sabnis,  et  al.  [6].  The  particle  trajectory  models  are  suitable 
for  simulation  of  steady  flow  problems  only,  while  the  particle  distribution  models  can  be 
used  for  steady  as  well  as  unsteady  flow  simulations.  The  algorithm  developed  by  Sabnis,  et 
al.  [6]  has  employed  transformation  of  coordinates  to  integrate  the  equation  of  motion  of 
the  particle  in  Eulerian  computational  space,  and  coupled  it  with  a  linearized  block  implicit 
scheme  of  Briley  and  McDonald  [7],  resulting  in  an  efficient  computer  code  for  two-phase 
flow  simulation  including  sub-layer  resolution  capability  without  the  need  to  use 
wall-functions. 

The  present  paper  describes  the  application  of  the  CELMINT  (Combined 
Eulerian-Lagrangian  Multidimensional  Implicit  Navier-Stokes  Time-Dependent)  code  of 
Sabnis,  et  al.  [6]  to  simulate  two-phase  flows  with  and  without  evaporation  and  comparison  of 
the  computed  results  with  experimental  data.  The  test  cases  used  in  data  comparison  were 
the  two-phase  mixing  layer  experiment  of  Milojevic,  et  al.  [8]  which  involved  no  evaporation, 
and  the  evaporating  spray  experiment  of  Solomon,  et  al.  [9]. 


REFERENCES 


1.  Gough,  P.S.:  Numerical  Analysis  of  a  Two-Phase  Flow  with  Explicit  Internal 
Boundaries,  IHCR  77-5,  Naval  Ordnance  Station,  Indian  Head,  MD,  April  1977. 

2.  Gibeling,  HJ.,  McDonald,  H.  and  Banks,  N.E.:  An  Implicit  Numerical  Analysis  for 
Two-Dimensional,  Two-Phase  Turbulent  Interior  Ballistic  Flows.  AIAA  Paper 
84-0561, 1983. 

3.  Crowe,  C.T.,  Sharma,  M.P.  and  Stock,  D.E.:  The  Particle-Source-in-Cell  Method  for 
Gas  Droplet  Flow.  ASME  Journal  of  Fluids  Engineering,  Vo!.  99,  No.  2, 

1977,  pp.  325-335. 

4.  Gosman,  A.D.  and  Ioannides,  E.:  Aspects  of  Computer  Simulation  of  Liquid- 
Fueled  Combustors.  Journal  of  Energy,  Vol.  7,  No.  6, 1983,  pp.  482-490. 

5.  Dukowicz,  J.K.:  A  Particle  Model  for  Liquid  Sprays.  Journal  of  Compu¬ 
tational  Physics,  Vol.  35,  No.  2,  April  1980,  pp.  229-253. 

6.  Sabnis,  J.S.,  Gibeling,  H.J.  and  McDonald,  H.:  A  Combined  Eulerian- 
Lagrangian  Analysis  for  Computation  of  Two-Phase  Flows.  AIAA  Paper 
87-1419, 1983. 

7.  Briley,  W.R.  and  McDonald,  H.:  Solution  of  Multidimensional  Compressible 
Navier-Stokes  Equations  by  a  Generalized  Implicit  Method.  Journal  of 
Computational  Physics,  Vol.  24,  No.  4, 1977,  pp.  372-397. 

8.  Milojevic,  D„  Borner,  Th.  and  Durst,  F.:  Prediction  of  Turbulent  Gas- 
Particle  Flows  Measured  in  Plane  Confined  Jet.  Proceedings  of  the  World 
Congress  on  Particle  Technology,  Nurnberg,  April  1986,  pp.  485-505. 

9.  Solomon,  A.S.P.,  Shuen,  J.-S.,  Zhang,  Q.-F.  and  Faeth,  G.M.:  A  Theoretical 
and  Experimental  Study  of  Turbulent  Evaporating  Sprays,  NASA  Contractor 
Report  174760,  September  1984. 

10.  Schlichting,  H.:  Boundarv-Laver  Theory.  Seventh  Edition,  McGraw  Hill  Book 
Company,  New  York,  1979,  pp.  729-750. 


THE  INTERIOR  BALLISTIC  FLOW  MODELLING 


AT  THE  ERNST- MACH- INSTITUTE 


R.  Heiser 


Fraunhof er-Institut  fUr  Kurzzeitdynamik 
-  Ernst-Mach-Institut  - 
Hauptstrasse  18,  D-7858  Weil  am  Rhein 


Abstract 


The  activities  in  interior  ballistics  modelling  at  the 
Ernst-Mach-Institut  are  presented.  Two  different  types 
of  flows  are  addressed: 

1 .  The  AMI  code  covers  the  multidimensional  two-phase 
flow  of  a  granular  charge.  Arbitrary  geometries  of 
the  flow  region  as  well  as  igniting  devices  on  the 
boundaries  are  included  now.  Eulerian  equations  are 
used . 

2.  The  multidimensional  turbulent  one-phase  flow  is 
investigated  by  an  extension  of  the  DELTA  code, 
an  Navier-Stokes  solver. 

For  both  cases,  the  theoretical  basis  as  well  as  some 
numerical  results  are  presented. 
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INVESTIGATION’  OF  UNSTEADY  FLOW  AND  COMBUSTION  PHENOMENA 


G.  Klingenberg 


ABSTRACT 

Prepressurizeb  mixtures  either  of  hydrogen/oxygen  or  hydrogen/ 
carbon  monoxide/oxygen,  with  helium  as  diluent,  are  filled  into 
a  windowed  combustion  chamber  of  a  single-stage  gas  gun  of 
caliber  20  mm  and  electrically  ignited.  The  interior  ballistics 
flow  and  combustion  phenomena  are  studied  by  means  of  pressure, 
temperature,  concentration,  and  gas  velocity  measurements. 
Infrared  emission  spectroscopy  is  employed  to  determine  local 
gas  temperature  and  concentration  profiles  of  chemical  species. 
A  multi-channel  spectrometer  is  used  to  measure  simultaneously 
the  IR  radiation  of  the  chemical  species.  The  results  presented 
provide  input  data  for  unsteady  fluid  dynamics  models. 


U.S.  Army  Workshop  at  Imperial  College 


Velocity,  pressure  and  heat  transfer  in  a  gas  gun 


R.  Perez-Ortiz 

The  evolution  of  gun  simulators  at  Imperial  College  has  led  to  the 
development  of  the  Mark  in  simulator  which  was  designed  to  increase  the 
projectile  and  gas  velocities  to  supersonic  speed.  The  Mark  III  rig 
simulates  a  constant-diameter  prebumt  propellant  ideal  gas  gun.  It  is 
designed  to  operates  horizontally  and  permit  initial  pressures  form  10  to 
60  bars. 


Gas  velocity  measurements  were  obtain  by  a  dual-beam  laser- 
Doppler  velocimeter  at  four  locations  from  the  breech  and  with  initial 
pressure  of  10  bar  gauge,  corresponding  to  muzzle  velocity  of  220  m/s. 
Indication  of  a  turbulent  wall  boundary  layer  thickness  ranging  from  26% 
to  43%  of  the  bore  was  found.  For  initial  pressure  of  10  to  60  bars 
corresponding  to  muzzle  velocities  of  220  to  420  m/s,  the  breech  pressure 
was  measured  by  a  piezo-electric  transducer  and  the  surface  temperature 
with  a  fast  response  thin-film  resistance  thermometer  at  seven  locations 
along  the  tube.  Comparison  of  measured  projectile  position  and  that 
calculated  by  assuming  quasi-steady  one  dimensional  isentropic  expansion 
shows  good  agreement.  Measured  temperature-time  profiles  at  different 
locations  and  the  same  pressure  show  that  the  temperature  change  due  to 
compression  heating  increases  with  distance  from  the  breech  and  the 
temperature  change  at  the  end  of  the  shot  reduces  at  locations  far  from  the 
breech.  At  the  same  position  the  wall-temperature  change  increases  for 
higher  pressures.  Preliminary  heat  flux  calculations  show  that  the  heat 
flux  is  high  immediately  after  the  projectile  passes  the  measuring  position 
and  subsequently  decreases  to  a  constant  value.  The  maximum  heat  flux 
occurs  at  locations  near  the  breech. 


Passive  and  Active  Combustion  Control 
K.  C.  Schadow  and  E.  Gutmark 
Research  Department 
Naval  Weapons  Center 
China  Lake,  CA  93555-6001 

Comoustion  characteristics  of  dump  combustors,  including  flammability  limits,  occurrence  of 
combustion  instabilities,  and  combustion  efficiency,  are  closely  related  to  the  shear-flow  dynamics 
associated  with  flow  separation  at  the  downstream-facing  step  into  the  combustor.  With  the 
detailed  understanding  of  the  shear-layer  development  from  initial  flow  instabilities  into  large-scale 
structures  and  subsequent  breakdown  into  fine-scale  turbulence,  it  is  possible  to  passively  and 
ac-jvely  control  the  combustion  characteristics.  Passive  combustion  control  is  achieved  by 
changing  the  initial  shear-flow  conditions  at  the  dump  using  nonstandard  inlet  duct  cross  sections. 
Active  combustion  control  is  obtained  by  acoustic  excitation  of  the  initial  shear  layer  and  periodic 
fuel  modulation.  The  control  experiments,  which  resulted  in  reduced  pressure  oscillation 
amplitudes  and  extended  flammability  limits,  were  conducted  in  three  steps  of  increasing  physical 
complexity.  Nonreacting  experiments  (hot-wire  anemometry,  flow  visualization)  were  followed 
by  combustion  experiments  in  laboratory  flames  (Planar  Laser  Induced  Fluorescence  (PLIF) 
imaging  technique)  and  dump  combustors.  The  combustion  control  experiments  with  dump- 
stabilized  flames  are  being  extended  to  disk-stabilized  flames.  Also,  in  addition  to  dump 
combustors,  the  critical  role  of  large-scale  structures  in  driving  pressure  oscillations  in  a 
regenerative  liquid  propellant  gun  is  being  investigated. 


PASSIVE  CONTROL  OF  COMBUSTION  OSCILLATIONS 

by 

Mario  Nery  R.  Nina 
Instituto  Superior  Tdcnico 
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An  experimental  investigation  has  been  conducted  in 
bluff-body  stabilized  premixed  propane  flames  in  pipes. 

Measurements  of  wall  pressure  were  performed  with  a 
Kistler  probe  located  in  the  cold  pipe. 

Results  obtained  with  different  pipe  lenghts  and 
stabilizer  geometries  are  presented,  bb  a  function  of  equivalent 
ratio  and  velocity,  in  order  to  characterise  the  domains  of  rough 
combustion  and  the  stability  limits.  The  results  obtained  with 
the  disc  were  used  as  reference  case. 
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ABSTRACT 


Intrinsic  burning  stability  of  solid  propellants  is  tightly  connected 
with  the  flame  structure  of  the  tested  propellant.  This  implies  the 
knowledge  of  the  fundamental  thermokinetics  not  only  in  the 
gas-phase,  but  also  at  the  burning  surface  and  often  in  the 
condensed-phase.  In  addition,  the  effects  of  the  energetic  sources 
driving  intrinsic  instability  may  be  significantly  increased  or 
damped  by  the  thermophysical  properties  of  the  propaoatmg  medium. 
Again,  this  implies  an  accurate  enough  knowledge  of  the  relevant 
properties  in  dependence  of  local  temperature  and  perhaps  pressure. 

stability  of  composite  and  double-base  solid  propellants  is 
discussed  within  a  nonlinear  framework  (finite  size  disturbances  of 
thermal  profiles  due  to  fluctuations  of  pressure  and/or  external 
radiation).  The  basic  assumptions  of  ID  and  quasi-steady  gas-phase 
are  invoked.  Specific  submodels  of  overall  thermokinetics  and  proper 
datasets  of  thermophysical  properties  are  presented  for  the 
condensed-phase  and  gas-phase  of  the  two  families  of  propellants 
under  etudy.  Double-base  compositions  are  less  stable  than  composite 
propellants  due  to  a  thick  fizz  zone  in  the  gas-phase,  extended 
reactions  in  the  condensed-phase,  and  also  less  heat  storage. 

Transient  burning,  computed  by  numerical  methods  under  a  variety  of 
operating  conditions,  allows  to  check  the  analytical  stability 
predictions  and  to  obtain  the  complete  time  solutions  of  the 
governing  set  of  equations.  Time-invariant,  self -sustained  steady 
oscillations,  and  the  trivial  unreacting  state  are  the  possible 
stationary  solutions  predicted  analytically  and  verified  numerically. 
Pressure  Deflagration  Limit  and  dynamic  extinction  are  also  predicted 
analytically  and  verified  numerically.  All  of  these  effects  are 
confirmed  experimentally  as  well.  Overstability  phenomena,  predicted 
analytically  and  verified  numerically,  were  experimentally  confirmed 
for  double-case  propellants  subjected  to  IR  laser  pulses. 

Combustion  wave  stability  is  augmented  by  heat  release  in  the 
gas-phase  and  diminished  by  heat  release  at  the  burning  surface 
and/or  condensed-phase.  Maximum  stability  is  obtained  if  energy  i° 
released  in  the  gas-phase  by  a  diffusion  flame  anchored  to  the 
burning  surface. 


A  Planar  Reacting  Shear  Layer  System  for 
Turbulent-Combustion  Interactions 

by 

C.J.  Marek,  C.T.  Chang,  B.  Ghorashi, 

C.C.  Wey,  C.  Wey,  E.J.  Mularz 

NASA  Lewis  Research  Center 

Several  recent  studies  and  reviews  of  current  research  in 
turbulent  reacting  flows  have  concluded  a  need  for  additional 
work  in  planar  free  shear  layer  systems.  In  the  mean  time,  there 
is  an  urgent  need  for  further  understanding  of  the  planar  mixing 
layer  flow-combustion  interaction  as  well  as  the  need  for  experi¬ 
mental  results  to  verify  the  accuracy  of  existing  computational 
fluid  dynamics  codes.  These  experimental  data  are  needed  for 
immediate  applications  in  space  and  aeronautics  projects.  For 
this  reason,  a  very  versatile  planar  reacting  shear  layer  (PRSL) 
facility  for  hydrogen  and  heated  air  (870K)  was  built  and  is  in 
early  testing  at  NASA  Lewis  Research  Center. 

The  multiple  objectives  are  (1)  to  establish  the  effects  of 
initial  conditions  on  jet  development  above  Mach  0.25  for  the 
prediction  of  the  shear  layer  spreading  rate,  (2)  to  identify  the 
effects  of  heat  release  on  the  turbulent  flow  field  of  the 
variable  density  flows  (i.e.,  density  difference  between  fuel  and 
coflowing  air),  and  (3)  to  investigate  the  effect  of  large  scale 
structures  and  acoustic  excitation  interaction  on  fuei-air  mixing 
in  high  speed  flow. 

Multiple  measurement  techniques  may  be  used  to  study  the 
flow  field  simultaneously.  Instrumentation  includes  laser  imag¬ 
ing  -  MIE  and  PLIF  to  determine  location  and  species  concentrat¬ 
ion  in  the  reaction  zone,  schlieren  -  time  resolved  density  field 
visualization,  LDV  -  flow  field  velocity  and  stress  tensor 
mapping,  and  acoustics  -  pressure  boundary  condition  data.  A 
description  of  the  facility  will  be  presented  along  with  the  data 
that  has  been  obtained  to  date. 


On  The  Flow  Structure  in  Two  Phase  Shock  Tubes 
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Abstract 

Numerical  methods  are  employed  to  investigate  the  structure  of  two  phase  Riemann  prob¬ 
lems  in  propellant  beds.  The  model  used  is  one  based  on  the  equal  pressure  model.  Since 
the  equations  arising  from  this  assumption  are  known  to  be  ill  posed  a  careful  analysis  of  the 
solutions  is  made  to  determine  how  this  lack  of  hyperbolicity  manifests  itself  in  the  results 

The  effects  of  particle  site,  gas  to  particle  density  ratio  are  also  investigated  and  the  results 
compared  to  both  dusty  gas  and  classical  shock  tube  flows. 

Finally  the  structures  calculated  are  examined  for  their  physical  plausibility  in  describing 
pressure  waves  and  deflagration  to  detonation  transitions  in  propellant  beds. 
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PROGRESS  IN  THE  CALCULATION  OF  UNSTEADY  AIRFOIL  FLOWS 
AT  LARGE  ANGLES  Of  INCIDENCE 


by 

Tuncer  Cebecl 
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California  State  University.  Long  Beach 


Abstract 


The  presentation  will  describe  the  progress  made  In  the  extension  of  the  steady 

Interactive  boundary-layer  method  of  Cebecl  et  al.  to  unsteady  flows  over 

practical  airfoils  subject  to  a  ramp-type  motion.  The  method  makes  use  of  the 
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unsteady  panel  method  developed  by  Platter  and  his  student,  Teng  .  and  Is 
able  to  compute  flows  with  large  regions  of  flow  separation.  By  solving  the 
quasi-steady  and  unsteady  boundary-layer  equations  In  an  Interactive  method, 
the  quasi-steady  method  will  be  assessed  over  a  range  of  angles  of  attack  and 
frequency  In  terms  of  convenience,  accuracy,  and  the  computational  cost-  The 
calculations  will  encompass  airfoil  and  wake  flows  at  angles  of  attack  close 
to  the  start  of  the  dynamic  stall  and  will  provide  Insight  into  the  development 
of  dynamic  stall  as  a  result  of  the  tral ling-edge  separation. 
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ABSTRACT 

Since  several  years  at  ISL  investigations  on  gun  tube  erosion  have 
been  carried  out  to  become  a  detailed  insight  in  the  action  of 
different  oropellar.t  charges  and  distinc  modes  of  additive  addi¬ 
tion  on  wear  and  erosion.  During  this  time  a  special  interest  was 
directed  towards  the  well-known  influence  of  so-called  addit'ves 
on  erosion  reduction.  In  this  paper  some  important  features  of 
recently  done  erosion  measurements  will  be  presented,  showing  the 
erosion  behavior  of  some  normal  and  LOVA  propellant  charges.  The 
LOVA  propellant  charges  investigated  have  been  receiveded  from  ICT 
containing  different  chemical  compositions.  The  effect  of  additive 
on  inhibiting  erosion  was  investigated  in  using  titanium  dioxide 
(T1O2)  or  wax  or  a  mixture  of  both  with  a  ratio  of  45%/55%  'Swe¬ 
dish  additive).  The  additive  was  placed  in  the  combustion  cham¬ 
ber  on  the  one  hand  as  a  "liner"  around  the  front  part  of  the 
propelling  charge  and  on  the  other  hand  in  a  "tablet  form"  gust 
behind  the  base  of  the  projectile. 

For  experimental  studies  of  gun  erosion  single-shot  experiments 
have  been  carried  out  in  a  test-gun-device  of  20-mm-cal iber .  For 
yielding  informations  about  the  eiosive  action,  surface  sensors 
have  been  used  for  measuring  the  surface  erosion  with  so-called 
Knoop  indentations  and  the  erosion  of  edges  or  cracks  with  cut  in 
circular  grooves.  Surface  analytical  and  metallographicai  investi¬ 
gations  of  the  surface  sensors  have  been  done  by  IFAM.  Besides, 
the  wall  temperature  was  measured  using  Ni-steel  thermocouples. 

The  erosion  measurements  have  shown  that  it  can  be  assumed  that 
one  parametei  influencing  the  erosive  action  of  a  propellant 
charge  is  its  heat  of  explosion.  It  ca  .ies  out  that  the  erosion 
increases  with  rising  heat  of  explosion.  The  erosion  and  the  in¬ 
ner  wall  temperatures  (depending  on  the  heat  of  explosion)  can  be 
related  to  each  other.  This  prediction  is  also  valid  for  the  LOVA 
propellants,  which  therefore  behave  with  regard  to  this  depen¬ 
dence  like  the  normal  propellants. 

Towards  erosion  reduction  the  best  results  have  been  gathered 
with  additive  insertion  as  a  liner  around  the  front  part  of  the 
propellant  charge.  Moreover,  the  liner  works  most  effective  when  a 
special  liner  material  (Ti02-support)  is  used  setting  free  very 
fast  the  additive  particles  during  the  interior  ballistic  cycle. 


Studies  on  secondary  combustion  of  gun  exhaust  plumes. 

H.  Mad) 


Abstract 

The  main  part  of  the  muzzle  flash  of  guns  is  due  to  secondary  combustion  cf 
reactive  propellant  gases  and  CO  with  ambient  air.  However  the  mechanisms 
which  are  tnought  to  be  responsible  for  this  reaction  as  well  for  inhibiting 
reactions  (e.g.  Carfagno's  model)  are  rather  hypothetical  and  not  really  veri¬ 
fied  by  experiments  on  real  guns.  The  present  paper  deals  with  the  conditions 
leading  to  the  ignition  of  secondary  flash  oc  a  7.62  mm  rifle  using  ammunition 
without  chemical  flash  suppressants.  In  particular  we  have  considered  the  beha¬ 
vior  of  propellant  gas  balls  included  such  with  mixed  air  along  their  pathways 
from  the  Mach  disk  through  the  intermediate  flash  zone  down  to  '•he  location 
where  ignition  occurs.  Along  these  lines  measurements  of  velocities  (by  LDA), 
temperatures  (by  spectroscopic  methods)  and  of  oxygen  contents  (by  use  of  an 
SF.-tracer)  were  performed. 
Additionally  temperatures  were  varied  by  using  gun  tuDes  of  different  lengths, 
oxygen  content  was  varied  by  shooting  into  a  closed  vessel  filled  with  air  of 
variable  oxygen  partial  pressure.  Ignition  was  recorded  by  high-speed  framing 
and  streak  photography.  Hence  ignition  times  could  be  provided. 

It  was  found  that  che  propellant  gas/air  mixing  along  the  flow  axis  was  rather 
gradually.  It  is  considerably  increased  by  the  presence  of  a  vortex  ring.  Tem¬ 
perature  and  oxygen  distribution  along  the  pathways  are  not  uniform.  There¬ 
fore  it  was  not  possible  to  arrive  at  an  analytical  expression  for  the  igni¬ 
tion-time  delay  in  terms  of  temperature  T  and  mole  fraction  a:  of  oxygen.  Such 
an  expression  is  well-known  for  the  C^/T^  reaction.  However  an  integral  form 
which  takes  into  account  the  time  histories  of  T  and  a  could  be  found.  It  in¬ 
cludes  a  parameter  which  is  equivalent  to  the  activation  energy  of  that  re¬ 
action  which  controls  the  rate  of  the  branching  chain  kinetics  of  the  secondary 
combustion.  The  experiment  yielded  an  activation  energy  which  is  considerably 
lower  than  the  activation  energy  of  the  ^2^2  f®30^00,  Th"s  we  conclude  that 
the  governing  elementary  reaction  for  the  ignition  of  the  secondary  flash  is 
different  from  the  elementary  reaction  H  +  0?  ■+  OH  +  0  which  at  a  hiqher  tem¬ 
perature  level  controls  the  H2/O2  combustion. ^Possible  other  reactions  will  be 
discussed  in  the  paper. 


